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SUMMARY 

The theory oj the steady-flow, linear type of crossed- 
field , d-c plasma accelerator, idealized for the ease 
of no wall friction and no transfer of heat as such 
to or from the plasma and specialized to the case of 
constant static temperature and constant cross-sec- 
tional area of flow, is developed from the it (dividual 
e ([nations of motion of the three components of the 
plasma . The results are more comprehensive , move 
nearly complete , and more accurate than previous 
results on this type of accelerator . The effect of the 
ion cyclotron angle ur, which is the product of the 
ion cyclotron frequency and the ion mean free time 
between collisions with neutral particles and which 
is proportional to the axial component of the ion slip 
velocity, on both Joule heating rate and accelerator 
length is included * n the results and is shown to be 
small only for values of about 1 ()~ 3 radian or less . 

INTRODUCTION 

The theory of the stead y-flow, linear type of 
crossed -field, d-c plasma accelerator, idealized for 
the case of no wall friction and no heat loss from 
the plasma and specialized to the case of constant 
static temperature and constant cross-sectional 
area of flow, is developed herein. The theory of 
the accelerator for the same case lias previously 
been treated (for example, in refs. 1 and 2) on 
the basis of the equations for the one-dimensional 
flow of the (‘enter of mass of the- plasma. The 
results so obtained have two faults. The first, 
which is lack of completeness, can be remedied by 
additional analysis based on the equations of 
motion of the electrons and the ions. For ex- 
ample, expressions for the axial component of the 


applied electric field, an important quantity that 
could not be obtained from the analyses based on 
the equation of motion of the center of mass, have 
been so derived and several are discussed and 
compared in reference 3. Furthermore, it is desir- 
able to have equations for certain additional in- 
formation, such as the magnitude of the deviation 
of the average ion velocity direction from the 
direction of the average neutral-particle velocity. 
Apparently, the equations for obtaining this and 
similar results have not been presented except in 
reference 2 and there the drag of the Coulomb 
forces was not taken into account. In the pre- 
vious analyses of the accelerator, a loss accurate 
form of the generalized Ohm’s law was used than 
is used in the present analysis. Because of the 
use of the more accurate form, the results of the 
present analysis are more nearly exact, but, as it 
turns out, the improvement in accuracy is sig- 
nificant only for sufficiently large values of the 
cyclotron angle a or of the ions. 

The present paper develops the theory of the 
accelerator from the individual equations of mo- 
tion of the three components of the plasma. Jt 
thus constitutes a single unified treatment that 
obtains the basic results that have been obtained 
previously but obtains them with somewhat 
greater accuracy, and also the additional results 
(such as axial component of the electric field and 
ion velocity vector) that come only from the 
equations of motion of the charged particles. It 
is a more comprehensive treatment than previous 
ones and therefore should serve to clarify certain 
aspects about which there may have been some 
uncertainly in the past. It makes rather easy the 
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taking into account of the long-range forces be- 
tween charged particles, does not require any 
consideration of (he components of the tensor form 
of the electrical conductivity, and shows the effect 
of ion slip on both Joule heating rate and acceler- 
ator length. 

In reference 2, for better insight into the plasma 
physics that is involved in the accelerator, an 
analysis from a microscopic, or particle, point of 
view was presented. A portion of that analysis 
is repeated herein to emphasize what the actual 
basic, physics of the acceleration process is and 
thus to provide a background for the analysis from 
the macroscopic point of view. 

SYMBOLS 

The rationalized mksq system of units is used 
herein. 


B 

c 2 
c p 

V 

E 

j 

k 

m 

M 

a 

V 

Q 

t 

T 

H 

V 

w 


y 

z 


magnetic induction, webers/nr 

constant, y A F u 2 -- r/r'r^r! ? dimensionless 
k I e‘ Iv 

specific lieat at constant pressure, jotiles/kg- 
°lv 


charge 4 on positive ion, <*ouIombs 

electric field strength, volts/m 

current density, amp/nr 

Boltzmann’s constant , joules/part iele-°K 

mass of part icle, kg 

Mach number, herein identical to M H 

number density, particles/m 3 

pressure, newtons/nr 

collision cross section, in 2 


time, set* 

tern jierat lire, °K 


dimensionless variable 


then 
1 ett 


velocity, m/sec 

velocity of ion in moving coordinate system, 


m/sec 

distance in axial direction, m 
distance in vertical direction, m 
collision rate per unit volume, m“ 3 sec -1 


a decree of ionization, — — — — - 

n-i + n i n 3 

j3 rate of flow of neutral particles per unit area, 
n^ tX) constant, particles/nr-sec 
7 ratio of specific heats 

e interaction parameter, newtons-m 2 -sec 


e 0 permittivity of vacuum, (3<iirX 10 9 ) -1 farad/m 
f dimensionless length, - r 

A" / 3 

A constant, — j dimensionless 

*23 

A quantity in expression for cross section 
for long-range encounters, dimensionless 

ju constant, a €l2 > dimensionless 

*21i 

p mass density, kg/m 3 

<r scalar electrical conductivity, mhos/m 

r mean free time, sec 

oj cyclotron frequency, sec -1 

Subscripts: 
e center of mass 

i final 

0 initial 

x component in ^-direction 

y component in //-direction 

i t j ndices indicat ing species 

1 elect ron 

2 'on 

J neutral particle 

Symbols defined as vector quantities by use of 
bold-face type are represented in italic, type to 
indici te scalar quantities. 

PRELIMINARY BACKGROUND MATERIAL 

A schematic diagram of a longitudinal section 
along a continuous-flow, linear type of crossed - 
field, d-e plasma- accelerator is shown in figure 1. 
The cathodes are shown in the upper wall and the 
anodes in tin 4 lower wall, the imposed magnetic 
field is directed normal to and out of the plane of 
tin 4 piper, and the electric field is in the plane of 
the paper and directed upward, perhaps at a 
slight angle to the vertical for reasons to be 
discussed subsequently. 


Cathodes--. 



Anodes-"'" 


Fkjure 1. -Schematic of plasma accelerator. 
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Concepts of the basic mechanism of the plasma- 
accelerator, that is, concepts as to the manner in 
which the charged particles in a three-component 
slightly ionized plasma accelerate the bulk of the 
plasma, can be developed from two points of view. 
From the microscopic point of view, the behavior 
of an electron and of an ion in electric and mag- 
netic fields that are normal to each other is 
examined. If there were no collisions, the paths 
of both ions and electrons would consist of a 
helical motion with a superposed drift. The two- 
dimensional projection of these paths in a plane 
normal to the magnetic field would be a seiies of 
cycloids that lie along equipotontial lines in a 
direction normal to both the electric and the mag- 
netic field. Because of the difference in mass, the 
size of the cycle differs for ions and elections, and 
because of the difference in the sign of their elec- 
tric charges, the sense of rotation is opposite for 
the two. But the average direction and magni- 


tude of velocity arc the same for the two, 


EXB 
B B ’ 


in the absence of collisions. There are, of course, 
collisions, and the effects and the results of these 
are discussed in the next section. 

From both the microscopic and t he macroscopic 
point of view, it will be seen that the accelerating 
force on the plasma derives principally from the 
ions which, between collisions, are accelerated by 
the fields and then, by collisions with neutral 
particles, drive the neutral particles along the 
channel. From the macroscopic point of view, 
the driving force is also shown to be the Lorenlz 
force on a current in a magnetic field, jXB per 
unit volume, where the current will be shown to 
consist principally of the flow of electrons. The 
microscopic method is used herein to provide 
insight into the physics of the process and the 
macroscopic method is used to provide quantita- 
tive results for use in designing accelerators. 


QUALITATIVE ANALYSIS FROM PARTICLE POINT 
OF VIEW 

Because of the large difference in the masses of 
electrons and ions, these two species behave 
diffeientlv in crossed fields in a three-component 
plasma; that is, they have different mean free 
times, different cyclotion frequencies, and different 
velocity directions and magnitudes. The desir- 
able behavior on the part of each is discussed first. 


Because of its comparatively small mass, an 
electron cannot impart much momentum at a 
collision with a neutral particle. It is therefore 
not undesirable that the electron make many 
cycles between collisions with neutral particles 
and thus have its velocity vector directed essen- 
tially at random just before collision. 

On the other hand, the ion should make, on the 
average, just a portion of a cycloid between colli- 
sions with neutral particles. It is perhaps well to 
reiterate here that the basic acceleration process 
is considered to be first the acquisition of addi- 
tional momentum from the electric field by the 
ion, then the transfer of this additional momentum 
to a neutral particle by collision, and then the 
equal distribution, on the average, of this addi- 
tional momentum over the vjn 2 neutral particles 
per ion. The portion of the cycloid that it is 
desirable for the ion to traverse is one for which, 
on the average, three conditions are satisfied. 

The first of those conditions is that, on the 
average, at collision the ion's velocity is a specified 
amount greater than the average forward velocity 
of the neutral particle. Satisfying this condition 
allows additional momentum to he imparted to 
the neutral particle at a specified rate. Part of 
this additional momentum goes into random 
motion of the particle and tends to raise the tem- 
perature of the gas. On the other hand, the cooling 
effect associated with acceleration of the plasma 
lends to lower the temperature. Thus, satisfying 
this condition results in a partial control over 
temperature. 

The second condition is that the collision occurs 
when, oil the average, the instantaneous direction 
of motion of the ion is in the direction of the axis 
of the channel. The reason for this condition is 
that it is desirable to drive the neutral particles 
axially along the channel and not toward a wall. 

The third condition is that the ion travels on 
the average along the axial direction — That is, in 
addition to moving axially at the time of collision, 
the average velocity of the ion should lie along the 
direction of the axis and not be directed toward a 
wall where on contact with an electrode the ion 
would be neutralized. 

Some of these three conditions are mutually 
contradictory, but usually not seriously so. 
First, one notes that the ions’ paths without 
collisions can be a series of prolate, common, or 
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curtate cycloids. A cycloid is generated by a 
point within, on, or without a rolling circle, the 
center of which moves with a constant velocity 


, , , 4 . EXB 

of translation - p - p » 

I> r> 


If the distance of the 


generating point from the center of the circle 


is less than the radius of the circle — the 

< (B B) 


point describes a prolate cycloid, as shown in 
figure 2. If t lie velocity of the generating point 
with respect to moving axes fixed in t he center 
of the circle is designated by w, then the instan- 
taneous velocity of the ion with respect to co- 
ordinates fixed in the laboratory is 


EXB 
B B 


w 


What is desired is that, on the average, the 
ion's path be the result of repeating a small 
section of this cycloid (fig. 2)- the section be- 
tween the two indicated points, a section for 
which the initial and the final directions are 
nearly the same and which approximates a straight 
line. This desired result is perhaps more easily 
seen in the hodograph plane, as in figure 3. The 


instantaneous velocity of the ion is 


EXB 
B B 


+ w. 


The neutral particle has the velocity v 3 before 


y 




collision. The ion has the initial velocity v 2t0 
in its free path after its previous collision and is 
accelerated to v 2t/ before the next collision. At 
collision, its velocity vector remains on the same 
small sphere on which it lay when its value was 
v 2 jr and the velocity of the neutral particle re- 
mains on the larger sphere. (These spheres have 
a common center and are represented in fig. 3 by 
circles ) These velocity vectors lie at the op- 
posite' ends of collinear radii of the two spheres, 
both 1 adore and after collision. If the scattering 
is isotropic in the center-of-mass system, then, 
on the average, the two velocities after collision 
lie at the mutual center of the spheres. Thus, 
on the average, at a collision the ion gives up 
the velocity increment it gained from the electric 
field b id ween collisions and returns to the init ial 
velocity in its free path, v 2(0 . By inspection 
of figure 3 it can be seen that, since the radii 
of the two spheres are in inverse ratio to the masses 
of the particles, the neutral particle on the aver- 


age ga ns 


m 7 


?n> 


m 3 


of the difference in velocity that 


exist ec before the collision. From figure 3 it can 
also hr seen that the velocity lost by the ion at a 

collision is the fraction — „ times the difference 
m 2 , 1 

m 3 + 2 


between its average velocity v 2 ^ which is equal 

to — V2 / ) il,u l the velocity of the neut ral 
particle v 3 . This result comes from use of the 
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hodograph and does not account for thermal 
motions. An unpublished analysis by Adolf 
Buscinann of the Langley Research Center has 
shown that, if thermal motions of the ion 
and the neutral particle are taken into account, 

the above fraction becomes approximately* — - — 

' — 2 ! +i 

or — — Thus, since the ion loses in velocity 
m 2 + m 3 

at a collision what it gained in velocity from the 
electric field between collisions, one can write 
that 


/ x m 2 m-i ( x j. 


Thus, the horizontal component of the electric 
field is 


E 


z' 


/ft 2 771 3 1 

m 2 T“ 7ft3 CT 23 


as was derived in reference 2. This result is re- 
peated here only for subsequent use in demon- 
strating that the expression for E x derived in a 
subsequent section agrees within a negligibly small 
amount with this result from the hodograph. 

The analysis of the accelerator on a microscopic 
or particle basis can be extended much further 
but to do so is not necessary as the analysis has 
been carried through on a continuum basis. It 
may be mentioned, however, that of the three 
desirable conditions discussed, the first two are 
adopted in the subsequent analysis. The first 
one is, in effect, adopted because constant static 
temperature is specified. Thus, the rate of accel- 
eration is specified to be such that the decrease 
in temperature due to acceleration is just com- 
pensated for by the increase in temperature due 
to Joule heating. The second condition is adopted 
in order that the driving force on the neutral par- 
ticles may he directed parallel to the axis of the 
channel. The third condition cannot then be met 
and the average ion velocity vector is directed at 
a slight angle above the axis, as can be seen from 
both figures 2 and J. Equations for determining 
the magnitude of this angle are derived in a 
subsequent section. 


QUANTITATIVE ANALYSIS FROM MACROSCOPIC 
POINT OF VIEW 

In order to analyze the plasma accelerator from 
a macroscopic point of view, the equations of 


motion, of energy balance, and of conservation of 
mass are written in terms of average velocities. 
The three equations of motion for the three species 
in a three-component plasma have been adapted 
from those given by Schlutcr in reference 4 and 
quoted also in reference 5 (section 48) with the 
addition of thermal-diffusion terms which, of 
course, are not needed for the present analysis. 
Equations of motion for the charged particles 
were derived in reference 2, but they do not in- 
clude the effects of Coulomb forces; therefore, 
Sehlii tor's equations, which do include these ef- 
fects, are used herein. Schltiter’s equations of 
motion, all other equations taken from references 
5 and fi, and all data have been converted to the 
rat ionalized mksq system of units. The equat ions 
of motion for the electrons, ions, and neutral 
particles are, respectively, 


/q//qVi (VVi) T- V P \ + w i ^3*13 (Vi y 2 ) T 1? (Yi y>) 
— — /tiKviXB) — (l) 

'%w 2 v 2 (vv 2 ) +VP2- i- /e>//3e>3(v 2 — v 3 ) -f ■w 2 w-,€ 2 i(v 2 — v,) 

= n 2 ( j ( v 2 X B ) + ihc E (2 ) 

W3M3Y3 (VV 3 ) +VP:i + nPhesi (v 3 — 

+ «3R-2€3:>(V 3 — V 2 )^() (3) 

Consider first equation (1). The two terms on 
the right-hand side express the force exerted on 
the electrons by the applied magnetic and electric 
fields. The last two terms on the left-hand side 
are the “friction drags” on the electrons caused 
by, respectively, encounters with neutral particles 
and encounters with ions. The drag is set pro- 
portional to the velocity difference and propor- 
tional to the number density of each of the two 
species. The proportionality factor is e ijt which 
thus becomes the “friction coefficient” or inter- 
action parameter between particles of species i 
and particles of species./. In a subsequent section 
is discussed the determination of the value of the 
e terms for (a) encounters of charged and neutral 
particles from kinetic* theory, quantum theory, 
and experimental values of mobilities of charged 
particles, and (b) encounters of charged particles 
from the theory of long-range encounters. 

The first two terms on the left-hand side of 
equation (1) will be neglected in the present anal- 
ysis. They are the term for the acceleration of 
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the electrons under the net effect of all the forces 
and the gradient of the partial pressure of the 
electrons. The neglect of these two terms can be 
justified by comparison of their magnitudes with 
the magnitudes of the other terms in the equation. 
A typical set of values, calculated from the results 
of the present analysis, shows that the order of 
magnitude of each of the last two terms on the 
left-hand side and of the two terms on the right- 
hand side is about 10 5 to 10 ° newtons per cubic 
meter, whereas the orders of magnitude of the 
acceleration term and the pressure gradient are 
respectively, 1 0 ~ 1 and 10 newtons per cubic meter. 

One further term that has been neglected is the 
time-dependent- variation, at a fixed location, of 
electron momentum, since only the steady-state 
case is considered herein. It may be observed 

that Sehluter's notation-^ indicates the substan- 
ce 

tial derivative^*? which is J*-? + These two 

quantities are neglected herein in the equations of 
motion of electrons and ions, but the second one 
is retained in the equation of motion of the neutral 
particles, where, instead of (v . v)v, t ho equivalent 
expression v(vv) is used. 

The equation of motion of the ions is very 
similar to that of the electrons. The same simpli- 
fications can be made. For a typical case, the 
acceleration term and the partial pressure gradient 
have magnitudes of about 10* and 10, respectively, 
and can be neglected in comparison with the other 
four terms, which are of the order of 1() 5 to 1(V 5 
newtons per cubic meter each. 

Equation (3) is the equation of motion of the 
neutral particles. Here again the term containing 

has been omitted because only the steady state 

is being considered. The acceleration and the 
pressure gradient terms, the first two terms on 
the left, are large enough to be retained because 
they are functions of // 3 and, since the present 
analysis is made to apply to a slightly ionized gas 
(whet her the ions are of the same original species 
as t he neutral part icles or are of some other species 
used as seeding material), n- s is much larger than 
w-i* Thus, in a typical case, the first term in 
equation (3) may be of the order of about It) 5 
newtons per cubic meter, the second term about 
IT, and t lie third and fourth terms about 10 5 to l(/\ 


As is pointed out in reference 5 (section 48) the 
cqun-tkn of motion for the center of mass of the 
plasma can be obtained by additg together the 
three individual equations of motion (oqs. (1) to 
(3)) wi h all terms retained in the equations. This 
addition yields, for /q — iu } 

py c ( w c )+v„=jxB 

In this addition, all the “friction” forces between 
various components of the plasma have dropped 
out, inasmuch as these forces cannot cause motion 
of the center of mass of the plasma — that is, c u 
must equal * Si . 

For ir ubsequent use, equations (1) and (2), with 
the first two (relatively small) terms in each 
neglected, are added. The result is 


n v «-3«i3 Oi — v 3 ) + n 2 n^(\, — v 3 ) 

= H x e (v> — Vi) x B = j X B 

for 

■*h=-n 2 (4) 

and wi h the conventional definition of j, 


j = /qr ( V 2 v,) 


0 r >) 


Previous analyses of the plasma accelerator used 
the equation of motion of the center of mass, to- 
gether with equations expressing conservation of 
mass and of energy of t he center of mass. A some- 
what d fferent procedure is followed herein, in that 
there are used, as the basic equations that describe 
the fluid dynamics of the problem, the equation of 
motion of the neutral particles, an equation ex- 
pressin : conservation of mass of the neutral parti- 
cles, and an energy balance equation that is written 
under lie assumption that all of the change in 
energy of the plasma occurs in the energy of the 
neutral particles. The self-suhsistent equations 
of motion of the charged particles are used some- 
what as auxiliary equations, both to permit solu- 
tion of the fluid-flow pro Idem and to obtain 
additional results. Only the ease of constant 
cross-si ctional area of the channel is considered; 
thus, c< nservation of mass of the neutral particles 
is expressed by 


or 


V-Kv a) = (> 

( 'on st ant = 


00 


l nder die assumption that any change in energy 
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occurs only in the energy of the neutral particles, 
the energy balance equation is 

"3«W V 3 ■ v T, + //;i/w 3 (v 3 • v 3 ) (VV 3 ) -E j 

The further restriction is made that the tempera- 
ture of the neutral particles is kept constant along 
the length of the channel. As was pointed out in 
reference 2, there are justifiable reasons for this 
restriction. At the present state of development 
of the art, it seems reasonable to expect that a cer- 
tain minimum temperature is required to main- 
tain a given degree of thermal ionization in the 
plasma. At this point it is also appropriate to 
mention that heat transfer to or from the plasma 
by conduction, convection, and radiation is 
neglected herein. The energy balance equation 
thus becomes 

n 3 m 3 (\ 3 * v 3 ) (vv a ) = E • j (7) 

under the condition that 

y 3= Constant (8) 

It may be mentioned that the restriction of con- 
stant static temperature simply means that the 
(‘fleets of Joule heating on temperature and of the 
Lorentz force on temperature are made to be 
equal and opposite. This restriction could be ap- 
plied in t he present analysis either by deriving an 
expression for grad 7 T 3 and determining the con- 
dition for which it vanishes, as was done in ref- 
erence 2, or by setting grad 7 3 equal to zero in 
the energy-balance equation and wherever else it 
appears. Both methods haul to exactly the same 
result. The latter method was used herein. The 
general gas law can be used, expressed as 

]h — 7 3 (9) 

When the degree of ionization a is small, if can 
be expressed as 

a—- -■■■=( 'onstunt (10) 

In order that the plasma will be accelerated along 
the axis of the channel and not toward a wall of 
the channel, the Lorentz force should be parallel 
to the direction of the axis. Thus, with the 2-axis 
of the coordinate system in the same direction as 

622979 - 02 2 


the channel axis, there results the condition 

v ltX ~v 2tA (11) 

The condition that the neutral particles flow in 
the direction of the 2-axis is that 

^.* — 0 (12) 

The Mach number of the flow of neutral particles 
is given by 

7 M * = /err (M) 

When both ./*- and //-components are included, 
equations (1) to (Kj) constitute 16 simultaneous 
equations in 16 dependent variables. The equa- 
tions with both x- and //-components are equations 
(1), (2), arid (5). (The //-component of eq. (J) 
is tantamount to eq. (11).) The 16 variables 
are the scalars n u n 3 , p 3y and T» the two com- 
ponents ouch of the vectors \ ly v 2y v 3 , E, and j 
and the 2-component of A/. (For simplicity, 

the vector notation and the subscript J are 
omitted herein from A/, and M is understood to 
mean M 3 .) In order that the number of variables 
may not exceed the number of equations, the 

quantity B is taken to be constant. Specifying 
that B rather than K y shall be constant appears 
to be advisable, if for no other reason than that 
the length of the accelerator varies approximately 
as the square of the final Mach .number, as is 
shown in reference 2; if K y were held constant 
and B varied in such a manner that the static 
temperature of the plasma remained constant, 
the length of the accelerator would vary with the 
fourth power of the final Mach number, as is 
shown in reference 7 . In a subsequent section, 
e 13 and e 23 will be shown to be, for given species, 
functions only of temperature and, since tem- 
perature is constant in the present analysis, so 
are these e terms. The quantity e l2 will be shown 
to be a function of the (constant) temperature 
and proportional to the logarithm of It 

can, therefore, to a sufficiently good approxima- 
tion, be considered a constant. 

The plan of the solution of equations (1) to 
(1 J) is first to find expressions for tile components 
of the velocities of the charged particles and for 
the components of the electric field, all in terms 
of v 3 (or A/). Then these results, together with 
the other equations, are used for eliminating from 
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equation (M) all dependent variables other than M 
The resulting equation is integrated to provide 
the variation of M (or v :i ) with .r. Thus, since 
tlu v other variables can be expressed as functions 
of v ;j , the variation with x of all the dependent 
variables is obtained. 

In order to derive expressions for the compo- 
nents of the velocities of the electrons and of the 
ions, equations (1) and (2), with the first two 
terms in each neglected, are used. The vector 
product of equation (1) and B is formed. From 
the resulting equation the quantity VjXB is elim- 
inated by means of (‘(juation ( 1 ) and v 2 XB is 
eliminated by means of equation (2). After col- 
lection of terms, there results the equal inn 

i(//,ei;i) 2 -t 2/^€i ;i /oe u > t c 2 IP\\ l 

4 ' K 1 3 1 2 d ' ^ 2 W »€ 1 2 ( € 23 “h € 1 3 ) J V 3 

— ^' 3 ^ 13 ^ ( v 3 X B) ^ 3 ^ 13 ^ E -j - f 2 ( Ex B) (14) 

Similarly, the vector product of equation (2) and 
B is formed. From the resulting equation the 
quantity v 2 XB is eliminated by means of eq na- 
tion (2) and ViXB is eliminated by means of 
equation (l). After collection of terms, the re- 
sulting equation is 

[ (^ 23 ) 2 4 2th Ntf 1 •>€•>;< - e' 1 /** ) Vo 

- / 1 1 U 3 € 12 ( € 23 ■ “ € 1 3 ) V !+[(?< 3 € 2 3) 2 
+ »l^3*12(«234“ €13)] V 3-b^3 c 2«<?E 

-r ( v 3 x B) + e* ( E X B ) ( 1 5) 

A new dimensionless variable u and new dimen- 
sionless constants are defined as follows: 

_«»e»3 

W “ ell 

x-X" 

€33 

n 1 €,0 

~ - 

*h e>3 


^ eli 


two equations results in 
| < w. 2 +2/m 2 + 1 ) (X 2 7< 2 4-2juX?/ 2 4 - 1 ) 4 m 2 '/ 4 < 1 — X) 2 ] v, 

I ( u 2 12 n n 2 + 1 ) (X 2 -T- fi-\ X/x ) 

fJLU 2 { 1 x)(l + /x~r A/x) j« 2 V;j 

— [ ( A - 1 - m "I • X>lx) 7 / 2 -f- X ] " (v 3 XB) 

~ | (X + m + Xyu ) w 2 i Xj E 

-H(l -(-m+Xm)« 2 + 1] (EXB) (10) 

| (u 1 + 2 /x w 2 + 1 ) (A 2 1/ 2 4 2X/x?r + 1 ) + /x 2 m 2 ( 1 — X) 2 ] Vo 
[ ( 1 4 m4" X/x) (X 2 ?/ 2 { 2Xg?/ 2 4- 1 ) 

( fi u * ( 1 — X) (X 2 + /x + \fj ) ] u 2 v 3 

+ |X(X+/x+Xm) ? F+ 1 ] (v 3 XB) 

4" f X ( X 4‘ m 4- X/x) ?/ 2 +1] E 

■ 1 ■ I (X 2 4- X/x 4" n) w/ 4 ■ 1 ] y^ L , ( E X B ) ( 1 / ) 

Fro 1 1 equations (lt>) and (17), the x- and //-com- 
ponents of Vi and v 2 can readily be shown to be 
given by the following equations: 

| [(ur+Zpu? 4- 1 ) (Xb/ 2 4 2\/x7/“4 1) 4V« 4 (1 — A) 2 ]r, * 

^ [ (tr + 2ji v? + 1 ) (X 2 + M 4 X/x) 

- /x7/ 2 (1 - X) (1 4 n \~\fx)\irr :iiX 

— 1 1 X+g+X/x) ?/ 2 4~A] 7^4 - [(14 'm4~X/x)?/ 2 4-1| ^ 

(IS) 

| ( U 2 + li M 7V 2 + 1 ) ( A 2 7/ 2 + 2 X/x U 2 + 1 ) + /X 2 w 4 ( 1 - X) 2 1 r 1 1 y 

“ [ (X 4" M 4‘ X/x) u 2 4~ X] ut\ x 
“I (X4- m 4X/x)?/ 2 4 X] yy /X 

-|(l-|- M -|-X/i)7/ 2 -l] (19) 

[ ( «* +S-MM* + 1 ) ■ +2X M w 2 + 1 ) +M 2 « ' ( 1 -X)*|r 2 ., 

- | ( 1 +M + M ( X 2 U. 2 +2X/UM 2 + 1 ) 

+m« 2 ( 1 -X) (X 2 +m+ Xm)]wV 3 .t 

+ I X ( X + fX 4" X/x) ?/ 2 4’ 1 J y£ />!r 


and are introduced into equations (14) and (1/5). 
Then simultaneous solution for Vj and v 2 of the 


+ 1 (A 2 4~ Xjtx 4 " m) m 2 "f" 1 ] y^ (20) 
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|<t/ a +2/m*+ \)(XHt 2 +2\^+ 1 ) 

= — [ A ( X + + X/x ) ?/ 2 + 1 ] w\ x 

+ [X(\+/*+X/i)-M a +l] K {J 

-[(X 2 +Xm+m)m j + 1|^ (21) 

Equations (18) to (21) are valid for general appli- 
cation when the accelerations and the pressure 
gradients of the charged particles are negligible. 
The next step is to obtain expressions for E x aid 
E l{ in terms of u and r :itX for use in eliminating E x 
and Ky from equations (18) to (21). The expres- 
sions for K x and K y are obtained herein for tin* 
special case of a plasma accelerator' in which the 
plasma temperature is kept constant and in which 
the Lorentz force is parallel to tire direction of the 
axis of the channel. First, an expression for E u 
in terms of r :iiZ (and M) is found in the following 
manner. Tin 1 equation just preceding equation 
(4) is used to introduce jXB into equation (8), 
equation (9) is used to eliminate V/>s in terms of 
Vth and V 7 a, equation (0) is used to replace Vh 6 
in terms of V*v 3 , equation (8) is used to eliminate 
vT-t, and equation (18) is used to introduce XI. 
The result is that equation (8) is transformed into 

y\f 2 — | 

W 3 W 3 (v 3 *v 3 ) - - !\f> - (V-v 3 )= jXB-Vy (22) 

By use of the condition (e(|. (11)) that the Lorentz 
force has no -{/-component, equations (22) and (7) 
become, respect ively, 

yAP — 1 <h\ r . n 

>h>n 3 i\x ,/j. —J»K O’O 

>h»)s r xx~ / ' (24) 


Then, by use of e<juation (2;')), 

y fV= ( 1 — X ) ulir 3x 

I (x + m + Xm) >i 2 + i | (7.1 /-’ — 1 ) 


( 20 ) 


It is of interest, to note tlmt anotlier sot of eondi- 
lions can also be used to obtain the same expres- 
sion for h\. The conditions are 


r — q 

J - ,r ~ay 

for which, from the {/-component, of (equation (8), 

Xrj , tf = ^ 2 , y 

r Phis relation , used in conjunction with equations 
(19) and (21), also results in equation (26). The 
two sets of conditions are, of course, physically 
equivalent, inasmuch as, if there is no Lorentz 
force in the {/-direction, then 4 also should be no 
neutral particle velocity in that direction and no 
change in /> 3 in that direction. 

When h\ and 1% are eliminated from equations 
(18) to (21) by tlie use of equations (20) and (25), 
the result ing equal ions for t he velocity components 
of the charged particles are 


1 


Xu 


V Wa.jr 


(27) 


A Ut\ x , 

r '^ - Xr <-" = [(x+ M +XM)tt*+ ij(ri/*_ 1 ) 

The last step in solving the equations is to 
integrate the ^-component of equation (8). If 
equation (6) is used in the first term, equations 
(9), (8), and (6) are used in the second term, and 
equat ion (11) is used in the third and fourth terms, 
equation (8) can be written as 


From ('quations (28) and (24) it follows that. 


K 


yXP 


yXP — l 


Hr*. 


(25) 


This relation was also derived in reference 2. 

To find an equation for E Xi (‘(plat ion (11) is 
again used. The left-hand sides of equations (18) 
and (20) are then equal. By equating the right- 
hand side's of tin* equations and effecting consider- 
able simplification, the desired relation is found to 
be 


(X+M+ Xju) W“ 1 
1 -X 


E x —u(E 9 -~ Hr A tX ) 


0m 3 


dx 


</* 


+ «J^(ei3+ e 23.) {f'3,x ^l,x) — fl 


If equations (10), (18), and (6) are used, then 


kT, 


«0(ei3 + € 23 ) 


(yd/ 2 — 1) 


”777“ _= 


(29) 


An expression is needed next for r li<r — r 3tJ in terms 
of r 3ir . Equations (27) and (2S) yield the needed 
relat ion: 


^ lr 


. ?, 3.ar 

I (X+m+Xm)m'+ 1](t-K 


! -l) 
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Thus, since by definition 
can l)e rewritten as 


A- 


«13 

— , 


€23 


equation (29) 


where 


«(1-, \)e- 


( 37 ) 


2oj96 a (l+X) , 
k'I\ 


| (\-i-ju+A/i) «•+ 1 \(yM~ — l) 2 
7.1/- 


<K 7-'/ 2 ) 


(30) 


Inasmuch as u is a variable, it must be expressed 
in terms of .1/. By defining a constant c' 1 as 


kT/-ir- 


(31) 


it can be shown t hat 



Also j" can l>c defined as 


, _2a/Wl +X) 
S_ k 7 a 


(32) 


(33) 


Then equation (30) becomes 


<n 


(y T/”*)* 


(/( 7 .I/ 2 ) (34) 


Integration yields 

= (X + ■ n -\ \n)c- jj7*I/ 2 — 2 !"g,(y.l/' 2 )—_ ^v,J 

+~ Kt*!/ 2 ) 2 — 47.I/ 2 -! 2 log, (7 -l/ 2 )) " 

Under the assumption that a: ----- 0 when yd/ 2 —!, 
this result becomes 


That ihe quantity a is the conductivity for the 
case of no magnetic field is shown in appendix A. 
Except for the last term on the right, equation (JG) 
is exactly the same result that is obtained in refer- 
ences 1 and 2, if in those references the constant 
of integration is so adjusted that x is zero for 
7*1 f 2 equal to unity. 

In older to arrive at the significance of the last 
term i i equation (dG), which does not appear in 
the results of previous analyses, the principal 
difference between the present and previous analy- 
ses must be considered. The main difference lies 
in the statement of the electrical power input to 
the accelerator that is used in the energy-balance 
(‘((nation. In the present analysis the electrical 
power input is E j as given by equation (7). In 
the previous analyses, the power input E j was 
assumed to be given by equation (A5! of appendix 
A, in vhich the last term is the rate at which work 
is done by the Loren tz force on t he center of mass 
of the plasma and the rate of Joule heating is taken 

j ^ 

to be — * Equation (An), however, is based on 

er 

equation (A4), an approximate form of Ohm’s 
law. To evaluate Joule heating more accurately 
requires some algebraic manipulation, but the 
results arc interesting. By equation (Al) of ap- 
pendix A, Ohm’s law is 

1 =E +! + X I^XBi+X^XB)] 

Scalar multiplication by j gives 


f (X-| m i X/i^Ty .!/ 2 — 2 log, ( 7-1 1 1 ) qjyd 


+5 [ (7-)/ ') " — 4 7 *1 / - 1 2 log, (7 -l/ 2 ) I 3 1 (35) 


If equation (35) is divided t I trough by (X ’f- 
lu T X m) c* and equations (31) to (33) are used, 
there results 


■la IP 

P Jp’S.r 


J- 7.1 1 1 — 2 log, (7.1/ 2 )- 


7.1 / 2 


+ 


1 


2(X+ m +Xm) 


r-Ur- 
)><■ L 


4+ 


2 log, (yM 2 ) 

7il / 2 


( 30 ) 


— = E • j - • ~ (j X B) • C v , + Xv 2 ) 

cr i ~r A 

By use of equations (5) and (J 1 ), 
h ‘ ~Kj,-j,Bv ux 

a 

- " foyjv Ju X jy O'l , X 1’2, x) 

- J Kjy —jy Ih\x— n x e B (y 2 > w - v , . v ) (r l<x - i\ r ) 

By use of equation (28), this equation becomes 
:l “ = = KJ,—j»Bvt.T-n,e li t- 1 — 

cr A 
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and by use of equation (27), 

^J»—j^.T= Jw --+nrn^t li (vx a -r^ i +e a (p i . t -P 3 , x y\ 

cr 

m 

This same result is derived in appendix B by a 
more physical line of reasoning. Equation (38) 
shows that the rate of Joule heating, which is the 
difference between the power input and the rate 
at which the Loren tz force does work on the neu- 
tral particles, is given by the two terms on the 
right-hand side. The first of these two terms is 
the conventional expression for Joule heating 
when there is no magnetic* field, and the quantity 
cr in this term is the scalar conductivity given by 
equation (37). The second term on the right- 
hand side can be easily understood from either of 
two slightly different points of view. First, it is 
recalled that the* condition has been set that the 
net current density j in the accelerator is in the 
//-direction. In the /-direction, the components 
of ion and electron velocities are equal, and thus, 
although there is no net current in that direction, 
there are actually two equal and opposite currents 
along that direction. If these currents are defined 
relative to the speed of the neutral particles, then, 
according to equation (-18), the Joule' heating due 
to the currents is 


or 


1 «iaO J i) rw+ 1 easO's)™* 
a a 


ae 2 


O'l) ref + 




O 2 ) r#-| 


where the subscript rel means relative to ^ 3iX . But 
from equation (37) for the ease that Coulomb 
forces, as accounted for by ju? can be neglected 
(since there is no relative velocity between ions 
and electrons in the /-direct ion). 


for electrons and 


ae' 1 


= 

*13 


<7 > — - 

€23 


for ions. Thus the Joule heating rate due to the 
/-components of velocity of the charged particles 

is --- for the electrons and for the ions. From 

O’ j (T 2 

a slightly different point of view, it is seen from 
equation (38) that the loss (to heat) due to the 


slip between the driving and the driven elements 
is a function of the velocity difference, just as it is, 
for example, in fluid-coupling power transmission 
devices. 

To return to a discussion of the last term in 
equation (30), it should be clear that this term 
originates from the heating due to the difference 
between the /-components of the velocities of the 
charged particles and of the neutral particles. 
Tims, if f J 2 .x— Li,* approaches 0, then the last term 
in equation (30) must also approach zero. That 
such is the case can lx* shown in the following 
manner. The last term in equation (30) is pro- 
portional to but it is seen from the equation 

immediately preceding equation (30) that, as u 
approaches 00 ? r, >z — r :itX approaches 0. 

There is another* way of stating the same result. 
Equation (41), which is derived in the following 
section on evaluation of the f terms, shows that, 

1 m.+ m-s 

— = 60-7 Toi 

u m 3 

Thus, the last term in equation (30) vanishes as 
co 2 t 2 3 approaches 0. The last term therefore is 
negligible when ojot 2 3 is small enough. “Small 
enough” turns out to be approximately 10 ' 
radian, as is shown in figure 4. For this figure, 
the last term in equation (30) was written, by 
using equation (32), as 

Three values of the coefficient of the quantity 
within the brackets were used in figure 4. One 
value, zero, corresponds to the less accurate 
results obtained in references 1 and 2, in which 
the additonal Jonh* heating term discussed in the 


2(\+^+x^)c 2 
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present paper was not taken into account, and 
thus corresponds to the case of o> 2 t 23 approaching 
0. Another value for the coefficient, 10 -3 , 
applies to the example conditions discussed 
subsequently herein, for which co 2 ro :r — 1.5 X It)' 3 
radians when the Mach number is 2. The third 
value for the coefficient, 10 2 , corresponds to 
co L »ro 3 = -5X 10" 3 radian when the Mach number is 2. 
As the Mach number increases, however, the 
density decreases and thus <a*Toz increases. At 
a Mach number of 12, for example, uoTta is six 
times as large as it is at a Mach number of 2. 
Figure 4 thus shows that the effect of the additional 
source of Joule heating on the dimensionless 
length of accelerator required to reach a given 
value of Mach number is large for the larger 
values of Mach number if aj 2 r L > :i is not small enough. 

A similar kind of result is obtained for the 
current density j y . From equation (28) it can 
easily he shown by use of equations (J2), (J7). 
and (41) that 

n x e (1 +\)ui\ T 

^ 1 (X I m I Aju) w 2 -) 1 K 7 J/ 2 — 1 ) 

__ vKt'x.r 

<T lit\ 

I { /// o - m? 

(X T ju ( Aju) /W 3 2 

m 

Here again, if is small enough, the equation 
reduces to that derived in reference 2. 

The angle by which the direction of the average 
ion velocity differs from the axial direction is of 

Vo 

interest. The tangent of the angle is -~ iV and 

V2 ,r 

tlius, from equations (27) and (28), is 
<\u = M/ 

1 t [<X ■ n \\(yM J -i) 

It is of interest at this point to verify two 
statements that were made previously. In view 
of equation (11), the Lorentz forc(‘ depends on 
the //-component of j. Furthermore, r ltJ/ can he 
shown to be much larger than v> , y . Thus, the 
equation immediately preceding equation (4) 
verifies the statement in the section “Preliminary 
Background Material” that the accelerating force 



on I ho plasma can be thought of either as being 
the Lcrentz force jXB, in which j is due prin- 
cipally to l lie flow of electrons, or as being due to 
impacts on neutral particles of charged particles 
that had been accelerated by tin* applied fields. 
Furthermore, use of equation (11) makes it easy 
to show, by comparing the .r-components of the 
last two terms on the left-hand side of equation 
(X) and by using the subsequently developed 
fact that € 13 is of the order of 1 percent of e 23 , 
that the accelerating force on the neutral particles 
due to electron impacts is only about 1 percent of 
that due to ion imparts. This result agrees with 
the description of the acceleration mechanism 
list'd in the microscopic analysis, that the ions 
produce most of tin* accelerating force on the 
neutral particles. 

It cun also be shown that one of the results of 
the present macroscopic analysis agrees quanti- 
tatively with the corresponding result obtained 
by use of the hodograph. From equations (27), 
(28), (2b), and (41) 

1 


,X 

[(X+jx+Xg) u 2 +l](Xil/ 2 — 1) 

Er 

(FA)?/// 

— ( fn> 2+fn 3 )o) 2 T^ Ex 
(1 — \)m 3 lt 

_ (m 2 + rn 3 ) r 23 e E r 
(1— A )m 3 m 2 

In tl e preeeding section, the hodograph repre- 
sentation for collisions was used to find that 




(m 2 +m 3 )T 23 eE z 

m 2 m A 


Tf the effect of collisions between electrons and 
neutral particle's is neglected, A vanishes and 
the two results agree exactly. 


EV ALUATION OF THE FRICTION FACTORS e u 

The riction factors that appear in equations 
(1) to ( \) are given by equation (hi. 4) of reference 
5 as 


4 Vtt m r \ tfij/ 


( 40 ) 
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It is worth noting that, since kinetic theory shows 
the binary collision rate per unit volume between 
species i and species j to be 


Z = 

^ l 7 *> 



m r | m 
m i ntj 



then 


nii-tnij ry 

'tl’i'lbj — *ij — /'if 

m i m j 


It should l)e noted that //,//;€ ,7, wliieh reference 
5 states lias the intuitive gas-kinetic meaning of 
the number of collisions per unit volume per unit 
time, must be multiplied by (m,+ to 

give that number. Inasmuch as the number of 
collisions per ion per unit time with neutral 
particles is 

m 2 + /«3 

72/3623 


the mean free time is 


To 3=7 


m 2 772 3 


(m 2 + 7/73)773623 

Since the cyclotron frequency of the ion is 

eB 


w 2 = . 


then 


m : 


m > e 


B 


OJoT 23 = 


m 2 i- m 3 7/3623 


ficulty is due to the Ramsauer effect, whi<*h causes 
a, large variation in Q VA near the lower energies. 
This variation is illustrated by figure 1.3 of refer- 
ence 8, for example, in which the collision prob- 
ability l\ is proportional to the elastic collision 
cross section Q vs of equation (40). The apparent 
difficulty in evaluating 613 lies in the rapid and large 
variation of Q x 3 with velocity. This difficulty can 
be circumvented by using experimental or theo- 
retical results 011 mobilities or drift, velocities of 
electrons. The method is simple. From equation 
(1), for vanishingly small values of electron accel- 
eration, electron pressure gradient , neutral-particle 
velocity, and magnetic field, the momentum bal- 
ance is 

3*13^1 — ( E 


Thus, 



(42) 


The quantity ~ is the electron mobility. Data ill 

the form of either mobility or drift velocity can be 
used for evaluating 613. 

For the friction factor between ions and neutral 
particles 623, no quantum (Ramsauer) effects arc 
important; nevertheless, since Q 2 3 may he a func- 
tion of velocity and data on Q,>z itself arc very 
scarce, [’('suits on mobilities, if they exist, can be 
used again to evaluate e> 3 : 


and, from the definition of u or equation (32), 


unr 23= 


7/2 2+ 272 3 U 


(41) 


This result was used in the preceding section. 

Equation (40) involves the binary collision rate 
as calculated from thermal velocities. Although 
the equation will be used herein to calculate only 
6j2 and the temperature correct ion to experimental 
data for e 23 , a check should be made to verify that 
the collision rate is determined by thermal veloc- 
ities and not the transport velocities of the 
particles; in other words, it must be shown that 
the thermal velocities are much greater than the 
relative transport velocities. This will be done 
subsequently for an example set of conditions. 

The numerical evaluation of the 6 terms poses 
some difficulties. There are three e terms to be 
evaluated. One of these is e 13 for the collision of 
electrons with neutral particles. An apparent dif- 


/ /«; 




tin *'■> 


(43) 


The other difficulty lies in evaluating e V2} the 
friction factor due to long-range forces between 
charged particles, because there appears to be no 
wholly acceptable theory by which the cross sec- 
tion Q ]2 can be evaluated. Reference o (section 
03) gives the following formula for Q x 2 which is 
based on SpitzeEs work (for example, ref. 0) and 
is simplified here for singly ionized ions: 


Q, 


(47re„/’7 T ) : 


log.. A 


(44) 


The difficulty is in the evaluation of A, which con- 
tains the limit at which tin* ( Ymlomb forces can 
be considered to act . Reference 0 uses the Debye 
shielding distance for which 


A — 1 27 r 
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Reference 5 (section 63), on the Imsis of some ex- 
perimental results on an arc, uses one-fourth the 
distance between ions and obtains 

A-— 4?r ( 4 *>) 

Numerical evaluation of the e terms for a typical 
set of conditions is desirable in order to establish 
orders of magnitude. A nitrogen plasma at 
4,000° K, with dissociation neglecte<l and seeded 
2 percent with cesium, is assumed as an example. 
The quantity e ]3 depends on the collision frequency, 
which depends on the relative velocity of the 
neutral particles and the electrons due to their 
thermal motions. On account of the large mass 
ratio between these two species, however, t lie- 
relative velocity due to thermal motions is de- 
termined almost wholly by the thermal velocity 
of the electrons. In equation (40), the reduced 
mass is essentially that of the electron. In this 
equation, therefore, T can be taken to mean T } 
and Q n can be considered to be a function of 7\. 
When €ia is evaluated from data by means of 
equation (42), tin 1 dependence of e 13 on 7 7 , can 
easily bo taken into account by using data ob- 
tained at the desired value of 7\, in this case 
4,000° K. As is well known, when electrons 
drift through a gas at a steady velocity in a con- 
stant electric field, as in the experiments that 
measure drift velocity, t lie temperature of the 
electrons (and here the term “temperature” is 
used to indicate the energy of agitational motion 
of tin 4 electrons, which are not in thermal equilib- 
rium with the neutral gas and which do not have a 
Maxwellian velocity distribution) exceeds the 
temperature of the gas by a factor that can be 
large even at moderate values of Kip. The data 
that are used in equation (42) are taken from 
figure 1 of reference 10. These data were ob- 
tained at a gas temperature of 293° K. The 
ratio of 4,000 to 293 is 13.6. Figure 7 of reference 
11 shows that the temperature of the electrons is 
13.0 times that of the natural gas at a value of 
Kip of 0.5 volt/(em)(mm Hg) or 0.37 volt-m/ 
newton. Figure 1 of reference 10 then shows 
that, at this value of Kip, the experimental value 
of r x is 5.1 X10 3 m/sec. However, the theoretical 
Xormand value shown in this figure may be 
preferable to t ho experimental value. In the 
experiments, the electrons at 4,000° K apparently 
excite rotation and vibration of the nitrogen 


molecules, which are- at a much lower temperature. 
(Evide ice for this process is presented and dis- 
cussed in ref. 9.) This process increases the drift 
velocity of the electrons but presumably would 
play a much smaller role in a plasma in which 
the electrons and the neutral gas were essentially 
in equilibrium at the same temperature of 4,000° 
K. Thus, the theoretical curve of drift velocity 
(curve labeled Norm and extrapolated to lower 
values of Kip) gives perhaps better values. (See 
fig. 1 of ref. 10.) This method yields a value of 
of about 3.3 X 10 a m/ser. Thus, 



1 0* new tons/' volt -see 


But 

*’l V[h^k J ;j 

K/p K 

t hus, 

~V i =2.2X 10 24 m-'-volt-'-sec" 1 

and 

e K 

e i r — - 7.3X10 41 n e w t o n - n i *-se< * 


The nobiity of cesium ions in nitrogen is given 
on pag * 407 of reference* 9 at 1S° C and 760 mm 
llgas 

^,=2.3r>Xl() -4 m 2 /volt-s(>c 

Thus, 

').9v'10 21 in '-volt '-sec' 1 

E 

and 

f K 

e 23 = - - v =2.7X 10 11 ncwton-m 2 -see 
n&i 


Tliis is the value of € £3 at 1N° C. The temperature 
variation of or of e 23 , or of in equation (40), 

has apparently not been determined for cesium 
ions in nitrogen either experimentally or theo- 
retieall v \ For the present use, therefore, the effect 
of temperature on thermal velocity only, and not 
on ^23, is taken into account and yields a value 
at 4,00 )° K of 

€23 — 1 .0 X 10” 40 newt on-m 2 -see 


(Fortunately, uncertainties in the value of 
have a negligible effect on the value of <r as calcu- 
lated from eq. (37).) 
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In order to calculate the friction factor e V2 for 
long-range encounters, A must be calculated. If 
the static pressure of the neutral particles is taken 
as 0.1 atmosphere, the temperature as 4,000° K, 
and the ionized seed substance as 2-percent mole 
fraction of the neutral substance, then, from 
equation (45), 

log* A =3. 03 

Actually, this value and the ion density and tem- 
perature from which it was calculated are slightly 
beyond the limiting conditions where the theory 
of reference 0 might be expected to break down. 
Equation (40) gives 

log, A -2.73 


the angle is 11°, and for a Mach number of 2 or 
more it is less than 1 minute. 

From equations (27) and (28) it can be shown 
that 

Vi . z — r :i , x = r 2 ' X —v 3 , * = 5 . 9 j ) i /sec 
= “080 in /sec 
v 2 , y — 0.50 m/sec 

The thermal velocity of the electrons 

r=( j =4.3X1 0 5 m/sec 

is enough higher than any of these other velocities 
to make equation (40) easily applicable. 


With the result from equation (45), the cross 
section for long-range encounters is, from equation 
(44), 

# I2 = (S.SX10" 17 m 2 
Then, from equation (40), 

€i2=3.2X 10“ 11 newton-m 2 -sec 
Thus, for the assumed example conditions 

X= -=7.3X10-* 

e 23 

fi—ct - -=6.4X 10“ 3 

e-i:\ 

From equation (37) the electrical conductivity is 


SUMMARY OF EQUATIONS FROM MACROSCOPIC 
ANALYSIS 

For convenience, the equations that have been 
derived as a result of the macroscopic analysis 
an* collected and listed in the present section. 

For t lie steady-flow plasma accelerator with 
constant cross-sectional area, constant static tem- 
perature, constant applied magnetic field, Lorentz 
force directed among the channel axis, small degree 
of ionization, no friction and no heat transfer, and 
7 XP greater than unity, the variation of A/ (the 
Mach number of the neutral particles) with dis- 
tance x along the channel (whore x is measured 
from the place where y XP is equal to unity) is 
given by equation (30) as 


cr — 720 mhos/m 

The dimensionless variable u, calculated from 
equation (32) for the example conditions, for B 
equal to 1 weber per square meter and a pressure 
of 0.1 atmosphere, is 




Tir 


- 1 15 


The quantity t o 2 t 23 for the same conditions is, 
from equation (41), 


2(7 B 1 


x=yM 2 2 log. 


1 


! ) — 


•XP 


2 (X+/i+X/i) 


JyXP 

us L 


. 21og f ( 7 A/ 2 ) , 3 1 
A- ^7 M 2 J 

=yM*-2 log, (tM 2 ) —^4 


I 

2(XH~/i+X/x)c 2 


[ (yM 2 ) 2 — 4 7 A P + 2 log, ( 7 A/ 2 ) +3] 


W2T23 = 1 . 5 X 1 0 ” 3 

The angle bet ween v 2 and v 3 , whose tangent is 
— f can be calculated under the assumption that 

Vo, x 

at the example conditions the Maeli number is 2 
and 7 is 1.4. Then it can be shown that- for^A/ 2 —! 


In equation (30) : 

(1) The quantity a is the electrical conductivity 
of the plasma when 7i=0 and r 3iX = 0 and is given 
by equation (37) as 

_ a* 2 (l+_X)_ 

(X + M + V) f 2 3 
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(2) The quantity u is given by equations (32) 
and (41) as 

1 eB 

u~~n^ n 

7tt«>+ m 3 

=* CO »To 3 

m 3 

and thus when cc 2 t 2 3 is small enough the last term 
in equation (36) can be neglected. 

8 2 € ^ 

(3) From equation (31) the quantity c 2 = ^ 2 ^ 3 2 * 

(4) The quantities e 2 3, X, and are “friction 
factors” or functions of them and can be evaluated 
from theory or experimental data. 

The variation of r 3tJ with x is found from equa- 
tion (36) and the following equations (eqs. (13), 
(12) t and (8)): 


The remaining variables are then given by equa- 
tions (25) to (28) and (39): 

yM 2 


E x 


K ^w-i lir ~ 

(\—\)uBi\ x 


[ (X -1 - ix + Xfx) u 2 + 1 ] (7 J f 2 — 1 ) 


, x — T ] x — " ^ r> , v ~t ^3 . x 


1-2,y — ” X?T f 


Xur 


3 . X 


X E x 


3v= 


v [ (X + ^ + Xfx) u 2j r 1 J (yAf 2 — 1) 1 — X B 

vBv 3fT 


_ 1+ (X+M+^)>r_ {yA/J V> 



» 3 ,„ = 0 

7 3= Const ant 

Then the number densities are obtained by use of 
equations (6), (10), and (4): 

w 3^3, z = & 

U\ = an 3 

and p 3 by use of equation (9): 

Pz — njc 1 3 

The quantity j x is liom equation (11) 

jx = 0 


CONCLUDING REMARKS 

A theoretical treatment of the steady-flow, 
linear, crosse d-field, direct-current plasma acceler- 
ator for inviscid, adiabatic, isothermal, constant- 
area flow has been developed from the equations 
of motion of the three components of the plasma. 
The results are idealized for the ease of no wall 
friction and no transfer of heat and specialized to 
the case of constant static temperature and con- 
stant cross-sectional area. The effect of the ion 
cyclotron angle to 2 7- 2 3 (where is the ion cyclotron 
frequency and r 2 3 the ion mean free time) on Joule 
heating rate and accelerator length is shown to be 
small ( nlv for values of about JO -3 radian or less. 

L A N ( ; LEV R K SE A HC H C E X TE R, 

National Aeronautics and Space Administration, 
La e; ley Air Force Base Va., June J } 1061 . 



APPENDIX A 


GENERALIZED OHM’S LAW 


The generalized form of Ohm's law is easily 
obtained from the equations of motion of the 
charged particles when the accelerations and the 
pressure gradients of the charged particles are 
neglected. Equation (1) is divided by ej 3 , equa- 
tion (2) is divided by e 23 , and the two resulting 
equations are subtracted and simplified. The 
result is 


i= E + -j- [(v lX B)+X(v lX B)] (Al) 

cr 1 -j- A 


where 

<xe~(\ -fX) 

(X+M+^m)ci>3 


(A2) 


The quantity <r is thus the scalar conductivity 
defined by Ohm's law for the case of no magnetic 
field: 

j = crE (A3) 

For cases where the accelerations and pressure 


gradients of the charged particles are negligibly 
small, equation (Al) is perhaps the most accurate 
of tlie generalized statements of Ohm's law in 
general usage. Less accurate statements are, 
however, frequently used. In reference 12, for 
example, the form that was used is exactly equiva- 
lent to equation (Al) with X set equal to zero in 
the right-hand side of the equation. Inasmuch as 
v 2 <v L and X<^<CL the approximation used in 
reference 12 is a very good one. In reference 2, a 
less accurate approximation was used. As Ohm's 
law, the equation 


J=E+v c XB (A 4) 

£T 

was used so that the power input to the accelerator 
was 


E-j = j j + ( j X B) - v c (AS) 

<J 
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APPENDIX B 


JOULE HEATING 


In the present appendix, the rate of Joule 
heating in the accelerator is derived by a some- 
what different- reasoning than that used in the 
body of this paper. 

The rate per unit volume at which the electric 
field does work on the electrons is the dot product 
of the force per unit volume exerted by the electric 
field on the electrons and the velocity of the elec- 
trons, that is, the dot product of equation (1) 
and Vi : 

u 1*3*13 (vj — v 3 ) -V! X •n-i%c !2 (v l — v 2 ) • Vi -= — n x e (E- v,) 

Likewise, the rate per unit volume at which the 
field does work on the ions is given by the dot 
product of equation (2) and v 2 : 

n 2 n ;i e ->3 ( v 2 — v 3 ) ♦ v 2 + n x 2 ( v 2 — v i ) • v 2 = n 2 e ( E • v 2 ) 

Adding these two equations yields 

n i n 3 c, :{ (V[ — v 3 ) • Vi X (v 2 — v 3 ) • v, 

X y*iw 2 c, 2 (v 2 — Vi) • (v,> — v,) = E j 

The next step is to obtain the rate at which work 
is done on the neutral particles. Addition of 
equations (I) and (2) gives an equation for jXB 
IS 


arid t he dot product of this equation with v 3 gives 
the rate per unit volume at which work is done 
on the neutral particles: 

th w 3 «i3 (Vi — v 3 ) • v 3 + w 1 n 3 e 2 3(v2 — v 3 ) • v 3 = j X B v 3 

Then, by subtraction, 

Ej-jXB •v 3 — /q/r-3 f 23(vi v 3 ) 2 

X Wi w«c ?a (v 2 — v 3 ) 2 + n x n 2 t i 2 (v 2 — vj 2 

By the use of equations (11), (12), (28), (5), and 
(37), this equation can be transformed to read 

^ yj y j\< z — dh ( 1 X X) e 2 3 { Vi, X C 3 , /) 2 

X X X) v i , w 2 
X ^1^2612(1 X X)V| t ,/ 

= n 1 n 3 ( 1 X X) 1 23 [ (r 2 , * — r - A , x ) 2 
X (XXm + Xm)'Y* 2 ] 

— / / 1 [ci 3 (r, tJ . C 3 , /) ** T C‘23 ( c 2 , r ? ’3 . z) 2 ] 



a 


which is the same as equation (38). 
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